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3       4 
The   H(p,y)  He  reaction  is  studied  for  incident 

proton  energies  from  8-30  MeV.   The  experimental  data  is' 

composed  of  excitation  functions  measured  at  the  angles  of 

55°,  90°  and  125°  in  the  center  of  mass  system  in  1  MeV 

steps  from  8-15  MeV  proton  energy  and  in  0.5  MeV  steps  from 

17-30  MeV  proton  energy.   These  data  measured  the  energy 

dependence  of  the  cross  section  and  its  deviation  from  the 

sin  8  angular  dependence  of  the  electric  dipole  (El) 

radiation. 

Additional  angular  distributions  were  taken  at 

selected  energies  in  the  range  from  13-30  MeV  with  up  to 

nine  angles  measured  between  20   and  14  2   in  the  laboratory. 


This  number  of  angles  was  taken  to  insure  statistical  signi- 
ficance for  determining  the  coefficients  of  the  Legendre 
expansion  of  the  differential  cross  section  through  a.. 
The  experimental  data  is  shown  to  be  consistent  with  the 
assumption  that  the  observed  radiation  is  primarily  El  with 
an  additional  electric  quadrupole  (E2)  component  and 
negligible  contributions  from  higher  order  multipoles.   The 
E2  radiation  is  shown  to  be  consistent  with  a  model  of  the 
E2  cross  section  composed  of  the  sum  of  a  direct  interaction 
component,  and  two  Breit-Wigner  resonances.   A  resonance  of 
width  12  MeV  at  a   He  energy  of  35  MeV  is  observed  in  agree- 
ment with  previous  experimenters  and  a  second  previously 

unobserved  resonance  is  observed  with  a  width  of  3.8  MeV 

4 
at  a   He  energy  of  39.5  MeV. 


INTRODUCTION 

Investigations  of  the  properties  of  the  4He  nucleus 
are  an  important  method  of  improving  the  understanding  of 
the  nuclear  force.   Consisting  of  four  nucleons,  -He  falls 
in  a  transition  region  between  describing  the  nuclear 
interaction  in  terms  of  nucleon-nucleon  interactions  and 
statistical  descriptions  of  many  nucleon  interactions.   The 
properties  of  the  He  nucleus  show  a  number  of  differences 
from  predictions  based  on  the  two-nucleon  interaction  and 
the  assumption  that  the  two-nucleon  forces  are  additive. 
Similar  difficulties  are  observed  with  models  that  are 
based  on  systems  with  a  large  number  of  nucleons. 

The  application  of  models  of  the  nuclear  inter- 
action to  4He  seeks  to  predict  the  observed  properties  of 

4 
this  nucleus.   One  of  these  properties  is  the  fact  that   He 

is  the  most  tightly  bound  of  all  the  nuclei.    Another 
property  to  be  considered  is  the  root-mean- square  (rms) 
charge  radius  measured  using  the  elastic  scattering  of 
electrons.2-4   Indications  of  the  presence  of  many  body 
forces  are  obtained  in  measurements  of  the  charge  distri- 
bution of  the  nucleus  by  electron  scattering.   The  charge 
distribution  deviates  significantly  from  a  Gaussian  shape.4 


Also  to  be  considered  in  a  description  of  He  is  evidence 
that  the  ground  state  is  not  a  pure  spin  zero  (S=0)  shell 
model  state. 

Many  models  have  been  presented  to  describe  the 
observed  properties  of   He.   Their  success  has  been  varied. 

In  general  the  models  proceed  by  incorporating  short  range 

fi  —  7 
correlations  and  non-central  forces  into  the  interaction. 

Further  information  on  the  importance  of  ground  state 

correlations  is  sought  from  the  angular  distributions  of 

a  •     9 

HHe   photonuclear  reactions. 

Several  measurements  of  the   He  photonuclear 
reactions  have  been  made.   A  number  of  inconsistencies  were 
observed  in  the  early  works,  some  of  which  have  been 
resolved  by  more  recent  measurements.   The  behavior  of  the 
angular  distributions  has  been  interpreted  to  give  suppor- 
ting evidence  for  the  presence  of  ground  state  correlations. 
In  addition,  evidence  exists  for  the  presence  of  resonant 

excited  states  in   He  whose  influence  should  be  considered 

13  14 

in  a  complete  description  of  the  angular  distributions.   ' 

Evidence  for  excited  states  in   He  is  presented  in 
the  analysis  of  many  experiments.   This  evidence  is  derived 
primarily  from  reactions  in  which  4He  is  formed  as  the 
compound  nucleus.     Evidence  for  several  states  has  been 
obtained  from  phase  shift  analyses  of  nucleon  scattering 
data.   Other  states  are  seen  in  inelastic  electron 

1  c     17 

scattering.   '     All  of  these  states  have  positive  energy. 


4 
There  exists  no  evidence  for  bound  states  of   He.   Isospin 

selection  rules  indicate  that  the  photodisintegration  of 

He  can  be  useful  in  the  investigation  of  negative  parity, 

j71  =  l",  isospin  T  =  1  states  and  positive  parity  J   =  2  , 

4         3 
T  =  0,  1  states.   Some  studies  of  the   He  (y,p)  H  reaction 

3       4 
have  exhibited    resonance-like  structures,  but   H  (p,y)  He 

measurements  have  failed  to  confirm  the  existence  of  the 

18—21  4         3 

proposed  levels.        Measurements  of  the   He  (-Y,n)  He 

TT      + 

reaction  present  evidence  for  a  J   =  2  ,  T  =  0  resonance  at 

22 

a  gamma  energy  of  E^  35  MeV.     It  will  be  shown  that  the 

present  work  provides  additional  evidence  for  the  existence 

4 
of  this  level  at  35  MeV  in   He. 

Most  of  the   He  photodisintegration  experimental 
efforts  for  excitation  energy  in  the  range  from  20-40  MeV 
have  been  devoted  to  the   He  (Y,n)  He  reaction.   The  early 
studies  of  the   He  (y,p)  H  reaction  suffered  from  diffi- 
culties due  to  poor  statistics.   For  energies  below  35  MeV 
where  the  triton  tract  becomes  too  short  to  measure,  the 
diffusion  chamber  and  cloud  chamber  measurements  have 

difficulty  distinguishing  protons  which  arise  from  multi- 

11  12 
particle  final  states.   '     Studies  of  the  inverse 

reaction  3H  (p,y)  He  have  been  conducted  for  limited 

excitation  energy  ranges  from  20  MeV  to  34  MeV.   In  the 

present  experiment  measurements  of  the  energy  and  angular 

dependence  of  the   H  (p,Y)  He  reaction  were  taken  for 

excitation  energies  from  26  to  42  MeV. 


The  proton  beam  from  the  Triangle  Universities 
Nuclear  Laboratory  Cyclo-Graaff  was  used  to  provide 
incident  protons  in  the  energy  range  from  8  to  30.5  MeV. 
The  beam  was  incident  on  a  5y  self-supporting  tritiated 
titanium  foil  target  containing  about  0.14  mg/cm2  of 
tritium.   The  emitted  photons  were  detected  with  a  thallium 
doped  sodium  iodide  Nal(TI)  detector  for  up  to  nine  labora- 
tory angles  between  20°  and  142°.   It  will  be  shown  that 
these  measurements  indicate  the  existence  of  a  previously 
unobserved  level  in   He  at  about  40  MeV  excitation  having 
a  J11   of  2+  and  a  width  of  about  3  MeV  in  the  center  of  mass, 

An  appendix  is  included  which  describes  a  peak 
fitting  program.   This  program  was  designed  to  extract  the 
low  energy  tail  of  the  Nal  detector  response  curve  from  , 
beneath  the  peaks  due  to  gamma  rays  from  transitions  to 
states  other  than  the  ground  state.   The  shape  of  the  Nal 

response  curve  was  based  on  the  measured  response  to  the 

3        4 
H  (p,y)  He  reaction. 


CHAPTER  I 
SUMMARY  OF  LITERATURE 

4 
The   He  nucleus  has  been  the  subject  of  intense 

interest.   The  resulting  body  of  literature  is  extensive 

in  both  experimental  and  theoretical  studies.   This  body 

of  literature  is  reviewed  in  a  number  of  recent  works, 

particularly  the  survey  by  Fiarman  and  Meyerhof  (1973)  .15 

Many  different  theoretical  models  have  met  with  qualitative 

success  in  describing  the  observed  ground  state  properties 

4 
of  He,  but  have  been  unable  to  satisfactorily  describe  all 

the  observed  properties.   Further  input  for  model  deter- 
mination is  provided  by  experimental  studies  of  the  excited 

4  4 

states  of   He.   While  there  are  no  bound  states  of   He, 

experimental  evidence  exists  for  virtual  states  as  shown  in 

Figure  1.1.   Particular  consideration  is  given  here  to  the 

photon  excitation  states  having  spin  and  parity  (J77)  and 

isopin  (T)  quantum  numbers  as  follows: 

jtt  =  1".,  t  =  1   and 

jtt  =  2  +  ,  T  =  0,1 
These  states  have  been  studied  by  photodisintegration  and 
the  inverse  photonuclear  reaction.   The  experimental 
results  to  date  have  not  definitively  established  the 


Figure  1.1    Excited  States  of   He.   The  energy  levels  of 
4He  and  the  general  form  of  the  cross  sections  for  the  four 
reactions,  which  form  He  as  the  compound  nucleus,  are  shown 
(Fiarman  and  Meyerhof   ) . 
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existence  of  any  of  these  states;  although  studies  of  the 

4He(Y,n)  He  reaction  indicate  the  existence  of  a  J77  =  2  , 

22 

T  =  0  state  for  an  excitation  energy  of   35  MeV.  *   In 

addition  to  the  question  of  the  existence  of  excited  states, 
there  are  fundamental  questions  regarding  the  observed 
shapes  of  the  photoproton  and  photoneutron  angular  distri- 
bution as  a  function  of  the  gamma  ray  energy.   The  (y,n) 
distribution  is  predicted  to  be  backward  peaked  while 
experimentally  it  is  observed  to  be  forward  peaked  for 

excitations  less  than  27  MeV  and  then  to  become  backward 

22 

peaked  between  27  and  35  MeV  excitation.     For  energies 

greater  than  35  MeV  excitation  the  angular  distribution  of 
the  (Y/n)  reaction  is  forward  peaked  and  the  asymmetry 
approaches  the  asymmetry  of  the  (y/P)  reaction  channel. 

The  (Y,n)  forward  peaking  may  be  interpreted  to  be  due  to 

23 

correlations  in  the  ground  state  v/ave  functions.  J 

The  experimental  observations  of  the  photodisinte- 
gration  reaction  have  limited  resolution.   These 
measurements  are  averaged  quantities  in  a  given  energy 
interval.  One  method  of  relating  these  averaged  quantities 
to  the  properties  of  the  Hamiltonian  and  the  ground  state 
wave  function  is  the  sum  rule  approach.24   The  theory  of 
sum  rules  has  undergone  extensive  development  starting  from 
the  study  of  atomic  photoabsorption  processes.   For  the 
electric  dipole  (El)  radiation  in  the  long  wave  length 


9 

approximation  the  integrated  cross  section  is  given  by  the 
sum  rule 


El 
El    =  /°°a(e)dE 
a 
int  (1.1) 

=  I      e*  <0|  [Zk,  [H,  Zk]]|0> 

1  60  NZ  (  =  60  MeV  -  mb  for   4He).24 


The  wave  vector  | 0 >  is  the  ground  state  wave  function,  H  is 
the  nuclear  Hamiltonian,  Zk  the  position  operator,  and  ek 
is  the  effective  charge  for  the  kth  nucleon.   The  equality 
in  the  last  term  is  satisfied  if  the  position  operator  com- 
mutes with  the  potential  energy  operator.   This  condition 
is  violated  by  Majorana  or  Heisenberg  exchange  forces  or 
velocity  dependent  forces.   Each  of  these  three  interactions 

would  tend  to  increase  o  .  ..   Electron  scattering  experi- 

mt 

ments  indicate  that  the  contribution  due  to  exchange 
currents  may  be  included  as  follows: 

El  25-27 

a  .    j.  =60  NZ(1  +  a),  a  =  0.40.  ~  (1.2) 

"A 

(=  84  MeV-mb  for  ^He) 

Under  the  assumption  that  the  ground  state  is  spatially 
symmetric  and  including  the  charge  radius  of  the  proton, 
the  Levinger,  Bethe  (28)  bremsstrahlung-weighted  cross 
section  is  related  to  the  nuclear  radius  by  the  following 
dipole  sum  rule: 


10 


El  r°°l      El  ,,..,_ 

°b      =/oE°        {L)dE 

=    0.096    ££=-    (R2    -    R2) 
A-l         z  p 

(    =    2.48   mb   for      He) 

where  R   and  R  are  the  nuclear  and  proton  rms  radii 

z      p 

respectively. 

F] 
The  quantity  a   (E)  referred  to  in  the  sum  rules 

is  the  total  cross  section  for  all  electric  dipole  photo- 

4 
disintegrations.  For  the   He  reactions  there  are  five 

competing  reactions.   The  yields  for  these  five  reactions 

were  measured  by  Gorbunov  and  co-workers.   '     A  170  MeV 

Bremsstrahlung  beam  incident  on  a  cloud  chamber  was  used 

in  the  study.   Their  results  for  the  relative  yield  of  the 

various  reactions  are  shown  in  Table  1.1. 

The  measured  angular  distributions  indicated  that 

the  El  interaction  comprises  %94%  of  the  total  integrated 

El 
cross  section  up  to  170  MeV.   Approximating  a    by  the 

total  cross  section,  integration  of  Gorbunov ' s  data  for 

the  reactions  of  Table  1.1  gives: 

aEl1   =  9  5  ±  7  MeV  mb  and  o    =  2.4  i  .15  mb 
int. 

This  value  is  larger  than  that  predicted  by  the  simple 

sum  rules  (Eq.  1.1-1.3)  indicating  contributions  due  to 

exchange  and  velocity  dependent  forces.   When  these  forces 

.  34 
are  considered  as  in  the  sum  rule  analysis  of  Quarati 
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for  the  shell  model  (Table  1.1),  the  experimental  cross 
section  up  to  170  MeV  excitation  exhausts  95%  of  the  sum 
rule. 

These  results  indicate  the  usefulness  of  sum  rules. 
This  usefulness  may  be  increased  by  the  removal  of  the 
approximation  that  the  total  cross  reaction  is  due  to  the 
El  interaction.   In  Table  1.2  are  listed  the  sum  rules 
commonly  appearing  in  the  literature  along  with  their 
representations  in  the  shell  model  formulism.   The  results 
of  these  equations  for  the  shell  model  parameters  are 

-i    _  i  -I  i 

S  =  1.75  fm2,  k-,.  =  -^-j-fm        =    36  MeV  and  k  „  =  ^fm 
El    5 . 4  E2    2 .  5 

4 
These  values  give  the  rms  charge  radius  of   He  to  be 

% 
<r.   >   =1.62  fm.   It  has  been  suggested  by  Leonardi  and 

38 
Lipparini    that  the  final  entry  in  Table  1.2,  a      ,  be 

-  3 

E2 

measured,  as  the  a     sum  rule  for  the  electric  quadrupole 

~3 

strength  is  dependent  only  on  the  ground  state  wave  function. 

These  features  of  the  four  nucleon  system  have 
received  extensive  theoretical  treatment.   Many  different 
models  have  been  successful  to  varying  degrees  in  descri- 
bing these  features.   The  common  factor  in  most  of  these 
analyses  is  an  attempt  to  develop  a  method  to  properly 
account  for  the  contributions  of  many  body  and  velocity 
dependent  forces  in  an  orderly  way.   Many  of  the  models 
are  based  on  "realistic"  potentials  that  are  derived  from 
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7  8  39  40 
the  nucleon-nucleon  scattering  phase  shifts.  '  '   '     These 

potentials  fit  the  on-shell  matrix  elements  equally  well. 
Their  differences  are  embodied  in  their  treatment  of  off- 
shell  matrix  elements.   One  such  potential  is  the  Sussex 

potential  which  was  used  by  Szydlik,  et  al.    to  calculate 

4 
the  spectrum  of   He. 

4 
Photodisintegration  studies  of  the   He  nucleus  span 

the  excitation  energy  range  from  threshold  to  500  MeV.   The 
reaction  is  characterized  by  the  asymmetries  in  the  proton 
and  neutron  angular  distributions.   At  low  excitation 
energies  the  asymmetries  exhibit  structure  which  has  been 
interpreted  as  being  due  to  interference  between  the  electric 
dipole  (El)  and  electric  quadrupole  (E2)  absorption  ampli- 
tudes.22.    For  excitations  above  100  MeV  the  (y,n)  and 

42 
(y,p)  asymmetries  approach  one  another.     The  absorption 

of  photons  of  energies  of  100  MeV  or  higher  is  expected  to 

occur  predominantly  by  the  quasideuteron  mechanism  discussed 

by  Levinger.    Recent  calculations  by  Gari  and  Hebach  account 

well  for  the  total  cross  section  and  asymmetry  in  the  range 

50  <  E  <  140  MeV.23 

4 
Many  measurements  have  been  made  of  the  various   He 

photodisintegration  processes.   An  overview  of  these  experi- 
ments is  shown  in  Table  1.3.   The  history  of  the  experimental 

41-44 
measurements  is  one  of  controversy.        The  early  measure- 
ments gave  the  general  character  of  the  energy  and 
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angular  dependence  of  the  excitation.   For  the  cloud  chamber 
and  diffusion  chamber  results  there  are  difficulties  relat- 
ing to  small  numbers  of  events  in  a  given  excitation  range 
and  in  analyzing  the  tracks  of  the  product  particles  for 
excitations  below  30  MeV.   Particular  difficulty  is  associ- 
ated  with  the  measurement  of  the   He  recoil  for  the  (y ,n) 
measurement.   All  the  photodisintegration  experiments  suffer 
from  difficulties  associated  with  the  spectral  distribution 
of  the  photon  beam.   These  early  measurements  indicated  the 
existence  of  a  number  of  structures  in  the  cross  sections, 
but  could  give  little  definitive  information  concerning  the 
(y,n)  assymmetry. 

The  measurements  by  the  different  experimenters  are 
in  general  agreement  concerning  the  magnitude  of  the  (Y/P) 
cross  section.   However,  the  (y,n)  cross  section  is  quite 
variable.   Charge  independence  considerations  support  the 
results  obtained  by  Gorbunov  that  the  (y,n) ,  (YfP)  cross 
sections  are  the  same.11,7  ~     The  measurements  are  in 

agreement  for  an  excitation  energy  greater  than  35  MeV. 

,  50-52     , 
The  results  reported  by  Berman,  et  al.        and 

Busso,  et  al.48,49,  give  a  lower  (y,n)  cross  section  for 
energies  between  23  and  35  MeV.   In  the  case  of  Berman' s 
results,  the  cross  section  would  agree  with  the  measure- 
ments of  Gorbunov    if  multiplied  by  a  factor  of  1.8. 
This  disagreement  has  been  clarified  by  a  series  of  experi- 
ments by  Irish,  et  al.53,54  which  show  that  the  variation  was 


22 


due  to  a  reduced  target  thickness  from  bubble  formation  in 

liquid   He.     Measurements  with  liquid  targets  by  Malcom, 

22 

et  al.     give  the  (y,n)  cross  section  in  agreement  with 

Gorbunov.   The  5  atm  diffusion  chamber  measurements  by  Busso 
agree  with  other  experimenters  for  excitations  above  30  MeV 
but  fall  below  the  (y,p)  cross  section  at  lower  excitation. 
Recent  1.5  atm  diffusion  chamber  measurements  by  Arkatov, 

C.C.  3 

et  al.     with  improved  accuracy  in  the  low  energy   He 

recoil  measurement  give  the  (y,n)  cross  section  in  agreement 

4 
with  the  (y,p)  cross  section.   These  results  for  the   He  (y,p) 

ind  (y,n)  cross  sections  give  experiment; 

charge  independence  of  the  nuclear  force. 


and  (y,n)  cross  sections  give  experimental  evidence  for  the 

79 


4 
The  existence  of  evidence  for  excited  states  of   He 

has  been  reported  by  many  experimenters  as  seen  in  Table  1.3, 
There  is  little  agreement  in  the  structure  seen  by  various 
early  experimenters.   The  recent  experiments  with  better 
statistics  exhibit  much  less  structure  than  seen  previously. 
The  experimental  cross  sections  are  shown  in  Figure  1.2.   The 
significant  features  are  the  peak  of  the  cross  section  for  an 
excitation  energy  of  26  MeV  and  a  broad  peak  for  28  MeV  exci- 
tation.  Peaks  in  this  region  are  consistent  with  the 
existence  of  two  states  with  angular  momentum,  parity  Jtt  =  1 

and  isospin  T  =  1.   These  states  correspond  to  a  mixture  of 

1  3    13 

the  spectroscopic  states  singlet   P,  and  triplet   P1- 


4  ... 

Figure  1.2     He  Photodisintegration  Cross  Sections. 

This  figure  is  a  compilation  of  the  results  reported  by 
the  indicated  authors. 
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Evidence  for  a  third  excited  state  at  35  MeV  excitation 

with  angular  momentum,  parity  J   =  2   and  isospin  T  =  0  is 

22 

obtained  from  (y,n)  asymmetry  measurements. 

The  (Y/n)  and  (y,p)  angular  distributions  are  given 
in  terms  of  the  Legendre  polynomial  expansion  as 

da  =  A   (1  +  Jak  Pk  {cos  6)).81  (1.4) 

dft    o      k  1 

The  coefficients  contain  contributions  from  photons  of 
multipolarity  L  for  k  <  2L,  k  even  and  k  <  L  +  L  ,  k  odd 
where  L   is  the  multipolarity  of  an  interfering  photon. 
The  equivalent  expansion  in  terms  of  sin  6  and  cos  9 
prevalent  in  the  literature  is 

4^  =  B(ct  +  sin2  6  +  6  sin2  0  cos  6+6  sin2  6  cos2  6  + 

dfi 

e  cos  8  +  •  •  • )  (1.5) 

The  coefficients  Ba,  B,  B£,  B6  and  Be  correspond  to  ampli- 
tudes and  cross  terms  for  multipolarities  of  0,  1,  and  2. 
The  respective  relations  between  the  coefficients  of  the  two 
expansions  are  shown  in  Table  1.4.   Electric  operators  in- 
volve only  space  coordinates  and  magnetic  operators  only 
spin  coordinates  in  the  central  force  approximation.   For  a 
pure  1S   ground  state  the  allowed  final  states  of  the  out- 
going particles  are   S,,   P.^  and   D2  for  angular  momentum 
changes  less  than  two.   The  admixture  of  D  states  in  the 
ground  state  indicated  by  variational  calculations  and 
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Table  1.4 

Relations  Between  the  Legendre  Polynomial  and 
sin  6,  cos  6  Angular  Distribution  Expansions 


B  =  Ao   (Zf  a2  "  I  a4  + 


a  =  1  +  a2  +  a4  + 


3       5 

2  a2    8  a4  +  *•• 


=  =|  a3  -  7a5  +  ... 
~2    a2  "  8  H4  +  '  ' 


jc   "35 

6  =  ~8  a4  +  •• 


-3       5  a   , 
~2  a2  "  8  a4  + 


e  =  a,  +  a,  +  ac  +  . . 
±     3  o 
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4         4  n 

He  (p,p)  He  measurements  complicates  the  interpretation  of 

the  coefficients  of  the  angular  distribution.   Possible 

additional  contributions  may  be  observed  as  an  increase  in 

the  isotropic  component  from  El  transitions  from  the  D  state 

to  P  and  F  states  and  a  contribution  from  the  cos  0  term 

which  arises  from  interference  between  El  and  non-spin-flip 

Ml  transitions.     The  angular  distribution  measurements  by 

Wait12   ,  Gorbunov 10    and  Meyerhof  13    yield  negligible 

values  for  the  coefficients  a  and  e.   Their  results  are 

dominated  by  the  El  term  B. 

The  asymmetry  (a  )  of  the  angular  distributions  about 

G  =  90°  can  be  defined  in  terms  of  the  difference  between  the 

cross  section  at  55°  and  125°.   This  quantity,  as  will  be 

shown  in  Chapter  II,  can  be  written  as 

P->  I  nn°\ 

3  (cos  55  )  co    „ 

a   =  a,  +  =: t C[-  'x  aQ  =  a   -  .68  a., 

s     IP,  (cos  55°)   3     1         3 

for  the  Legendre  polynomial  expansion  (or  equivalently  the 

term  3  for  the  sin  6,  cos  6  expansion)  where  P2  (cos  55  )  =  0 

eliminating  a„.   These  terms  arise  from  interference  between 

the  electric  dipole  and  quadrupole  interactions.   At  low 

energy  the  observed  angular  distributions  are  very  close  to 

a  pure  sin  6  distribution. D    This  distribution  indicates  a 

predominately  dipole  interaction.   The  isotropic  terms  were 

observed  by  Gorbunov  to  be  small  for  photon  energies  up  to 

170  MeV  providing  a  basis  for  the  assumption  that  only  El,  E2 
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radiation  contribute  significantly  to  the  yield.   The 
interference  terms  producing  the  asymmetry  are  more 
sensitive  to  smaller  amounts  of  E2  interaction  than  the 
E2  amplitude  itself  under  these  conditions.   This  result 
is  due  to  the  fact  that  the  coefficients  a,  and  a_  are 
proportional  to  the  product  of  the  dipole  and  quadrupole 
amplitudes . 


CHAPTER  II 
EXPERIMENTAL  METHODS  AND  MEASUREMENTS 

This  study  consists  of  a  series  of  experimental 
measurements  designed  to  obtain  the  electric  quadrupole  E2 
strength  as  a  function  of  the  excitation  of  the  compound 
nucleus  for  the  reaction  3H(p,Y)4He.   The  methods  used  to 
determine  the  E2  strength  consisted  of  measuring  the  asym- 
metry of  the  angular  distribution  about  90°  and  the  measure- 
ment of  the  detailed  angular  distributions  of  the  outgoing 
T-rays.   Shell  model  calculations  predict  a  number  of  states 
of  angular  momentum  and  parity  J11  =  2   with  isospin  T  =  0,1 
in  the  energy  range  from  30-50  MeV  excitation.   The  proton 
energies  available  from  the  Triangle  Universities  Nuclear 
Laboratory  (TUNL)  Cyclo-Graaff  allow  study  of  the  energy 
range  from  22  to  42  MeV  excitation  of  the  compound  nucleus. 
Experimental  Details 

Using  a  15  MeV  isochronous  Cyclotron  as  a  source, 
the  TUNL  FN  Tandem  provides  a  pulsed  25  MHz  proton  beam. 
The  beam  energy  is  continuously  variable  from  17-31  MeV. 
The  beam  pulse  width  is  2  ns  and  the  energy  spread  is  30  KeV. 
The  root-mean-square  beam  current  delivered  to  the  target 
was  kept  in  the  range  of  20  to  100  nanoamperes  depending  on 
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30 


the  total  count  rate  in  the  sodium  iodide  detector.   The 
count  rate  in  the  detector  was  kept  below  350k  counts  per 
second  because  pile-up  problems  and  a  deterioration  of  the 
resolution  occurred  for  higher  count  rates. 

Two  separate  detector  assemblies  were  used  in  the 
experiment.   The  first  detector  assembly  was  a  high  resolu- 
tion system  and  was  used  for  the  yield  curve  measurements 
and  the  angles  of  the  angular  distributions  between  42   and 
142°.   The  second  detector  assembly  allowed  measurements  to 
be  obtained  for  angles  between  10°  and  90   although  with 
reduced  resolution.   The  experimental  arrangement  was 
designed  in  both  cases  to  minimize  the  background  events 
that  were  not  associated  with  the  target. 

The  detector  assembly  is  based  on  a  cylindrical  . 
25.4  x  25.4  cm  thallium  doped  sodium  iodide  crystal.   This 
crystal  is  surrounded  by  a  well  type  Ne  110  plastic  scintil- 
lator.  The  center  crystal  is  viewed  by  six  RCA- 857 5  photo- 
multiplier  tubes  and  the  plastic  scintillator  by  eight  XP1031 
photomultiplier  tubes.  *"   A  schematic  of  the  detector  arrange- 
ment is  shown  in  Figure  2.1.   Active  shielding  is  provided 
by  detecting  coincidences  between  the  plastic  detector  and 
the  center-crystal  detector  which  were  processed  by  the 
electronics  to  reject  cosmic  ray  events.   Passive  shielding 
for  the  detector  is  provided  by  four  inches  of  lead  and  eight 
inches  of  paraffin  loaded  with  lithium  carbonate  (^50%  by 
weight) .   The  detector  is  mounted  on  a  carriage  that  allows 


Figure  2.1   Detector  Assembly  for  25.4  x  25.4  cm  Cylindrical 
Nal(Tl)  High  Resolution  Detector.   This  illustration  shows 
the  arrangement  of  the  active  and  passive  shielding  from  the 
cosmic  rays  and  beam  associated  background.   The  collimator 
arrangement  that  was  used  to  restrict  the  cone  of  gamma  rays 
to  that  just  illuminating  the  back  face  of  the  detector 
crystal  is  shown  (Weller,  et  al.82). 
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radial  and  angular  positioning  about  the  target  location. 
Two  radial  positions  of  the  detector  were  used  in  the 
experimental  series.   These  distances  placed  the  rear  face 
of  the  Nal  crystal  at  the  positions  of  82  cm  and  106  cm 
respectively.   Collimators  were  provided  which  restricted 

the  observed  solid  angles  to  those  subtended  by  the  rear 

-2 
crystal  face  corresponding  to  solid  angles  of  7.42  x  10    Sr 

-2 
and  4.46  x  10    Sr  respectively.   Excitation  curves  were 

measured  with  the  detector  at  the  82  cm  distance  in  the 
interest  of  higher  yield  rates.   The  angular  distributions 
were  measured  with  the  detector  at  the  106  position  to  allow 
the  inclusion  of  measurements  of  the  angles  at  42   and  142  . 
To  determine  the  angular  distribution  as  the  labora- 
tory angle  approaches  zero,  additional  measurements  were, 
taken  for  the  laboratory  angle  of  20   using  a  separate 
detector  system.   This  system  was  based  on  a  10.2  x  17.8  cm 
Nal(Tl)  crystal  surrounded  by  a  22.9  x  22.9  cm  Nal(Tl) 
annulus  for  active  cosmic  ray  rejection.   The  experimental 
arrangement  is  shown  in  Figure  2.2.   Passive  shielding  is 
provided  by  a  20.3  cm  thick  copper  annulus  and  4.4  cm  of 
lead.   A  cadmium  sheet  was  located  between  the  lead  and 
copper  for  thermal  neutron  absorption  and  the  collimator  was 
filled  with  paraffin  loaded  with  lithium  carbonate.   The 
distance  from  the  target  to  the  crystal  front  face  was 
93.3  cm.   The  collimator  restricts  the  photon  flux  to  the 


Figure  2.2   Detector  Assembly  for  10.2  x  17.8  cm  Cylindrical 
Nal(Tl)  Detector  for  Small  Angle  Measurement.   The  active 
and  passive  shielding  for  cosmic  rays  and  beam  associated 
background  is  shown.   The  beam  stop  shielding  for  measurements 


at  8,  ,  =  20°  is  shown, 
lab 
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6.57  x  10    Sr  solid  angle  subtended  by  the  back  face  of 
the  crystal.   The  entrance  aperture  was  formed  from  tungsten 
for  attenuation  of  fast  neutrons.   The  remainder  of  the 
collimator  was  lined  with  lead. 

Further  measures  to  minimize  the  background  events 
not  associated  with  the  target  consisted  of  the  use  of  a 
single  insulated  collimator  close  to  the  target  and  adequate 
shielding  of  the  beam  dump.   The  insulated  collimator  allows 
monitoring  of  the  beam  current  striking  it.   This  current 
was  kept  to  a  minimum  by  adjustments  of  the  focusing  of  the 
beam  transport  system.   The  beam  stop  for  the  large  crystal 
arrangement  is  located  three  meters  beyond  the  target. 
Originally  this  beam  stop  was  shielded  by  about  10  cm  of 
lead  and  70  cm  of  paraffin.   The  shielding  was  changed  to 
about  10  cm  of  copper  and  4  0  cm  of  lithium  loaded  paraffin 
to  more  effectively  attenuate  neutrons  produced  in  the  beam 
stop.   For  the  small  crystal  assembly,  the  beam  stop  was 
located  33  cm  behind  the  target.   The  shielding  was  provided 
by  the  copper  shield  of  the  detector  itself  and  a  tungsten 
shadow  bar. 

The  electronics  for  the  two  detector  systems  were 
similar;  however,  the  use  of  the  Nal  annular  cosmic  ray 
detector  on  the  smaller  crystal  required  a  different  elec- 
tronic arrangement  than  the  plastic  shield  of  the  25.4  x 
25.4  cm  system.   The  electronic  arrangement  for  the  large 
crystal  is  shown  in  Figure  2.3  and  that  for  the  small  crystal 


Figure  2.3    Electronics  Block  Diagram  for  High  Resolution 
Detector.   This  diagram  illustrates  the  signal  routing  used 
to  determine  the  coincidences  between  the  active  shield 
detector  and  the  center  detector.   The  additional  electronics 
used  to  determine  the  number  of  coincidences  are  not  shown 

(Weller,  et  al.82) . 
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in  Figure  2.4.   Typically,  the  count  rate  due  to  low  energy 
events  is  several  orders  of  magnitude  greater  than  the  count 
rate  due  to  events  of  interest.   Common  to  both  systems  is 
the  arrangement  of  the  electronics  to  minimize  the  events 
due  to  the  pileup  of  lower  energy  events.   This  minimization 
is  accomplished  by  time-clipping  the  detector  signals.   The 
signal  is  then  passed  through  a  fast  linear  gate  (250  ns) 
before  processing.   During  the  10  ys  processing  time,  the 
gate  is  held  closed. 

Further  reduction  of  the  background  due  to  unwanted 
events  was  possible  in  this  experimental  series  due  to  the 
pulsed  nature  of  the  Cyclo-Graaff  beam.   A  time  of  flight 
requirement  (TOF)  was  set  up  between  a  signal  derived  from 
a  radio-frequency  inductive  pick-up  loop  at  the  Cyclotron 
and  the  gamma  ray  events  in  the  center  crystal.   The  R-F 
signal  was  used  for  the  stop  signal  for  a  time  to  amplitude 
converter  (TAC) .   The  TAC  start  signal  is  derived  from  the 
linear  signal  of  the  Nal  center  crystal.   The  TAC  was  gated 
by  the  pile-up  rejection  circuitry  to  be  active  only  for 
the  time  during  which  the  linear  gate  was  open.   Both 
detectors  are  pulse  stabilized  by  a  calibrated  light  emit- 
ting diode  flash  into  the  crystal.   This  pulse  was  used  to 
correct  for  slow  changes  in  the  system  gain  due  to  thermal 
drifts.   The  pulse  stabilization  was  found  to  be  unnecessary 
for  the  large  detector  system.   The  cosmic  ray  rejection  was 
accomplished  differently  for  the  two  detector  systems.   The 


Figure  2.4  Electronics  Block  Diagram  for  Small  Angle 

Experiment.  Illustrated  here  are  the  electronics  used 

to  determine  the  coincidences  between  the  active  shield 

detector  and  the  center  detector  for  the  yield  measurement 

at  6.  ,  =  20  .   The  electronics  associated  with  the  pulse 

stabilization,  accidental  monitoring  and  beam  time  of 
flight  circuits  are  also  shown. 
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response  of  the  plastic  scintillator  is  faster  than  that  of 
Nal.   For  the  large  detector  with  the  plastic  shield,  the 
shield  signals  were  delay  line  clipped  and  timed  for  an 
overlap  coincidence  with  the  center  crystal  events.   When 
a  coincidence  occurs,  a  routing  signal  is  sent  to  the 
computer  for  storage  of  the  event  in  a  coincidence  spectrum. 
The  small  crystal  system  utilizes  the  linear  signal  from  the 
center  crystal  as  the  start  signal  for  a  time  to  amplitude 
converter.   The  shield  pulses  are  then  the  stop  signal  for 
the  TAC.   The  TAC  output  goes  to  two  single  channel  analyzers, 
one  of  which  is  set  to  provide  a  signal  for  events  occurring 
in  the  peak  of  the  coincidence  TAC  spectrum.   The  other 
analyzer  is  set  off  the  peak  to  determine  the  counts  due  to 
accidental  coincidences.   In  each  case  the  resulting  signals 
are  processed  by  two  analog  to  digital  converter's  (ADC's) 
in  coincidence  and  passed  to  the  on-line  computer.   For 
maximum  resolution,  the  threshold  of  the  coincidence  circuit 
could  be  set  to  give  a  coincidence  for  events  due  to  annihi- 
lation y-rays  from  events  in  the  center  crystal. 

The  information  stored  in  the  computer  consisted  of 
eight  512  channel  spectra.   A  window  light-penned  on  the 
y-flash  in  the  TAC  spectrum  separated  the  Y-ray  events  into 
three  pairs  of  spectra.   These  pairs  were  stored  as  short 
time  events,  true  data  events  and  long  time  events.   Each 
pair  corresponds  to  center  crystal  events  and  coincidences 
respectively.   The  eighth  spectrum  stored  was  used  with  the 
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large  crystal  system  to  store  the  elastic  proton  spectrum 
from  a  500y  silicon  surface  barrier  detector  located  at  a 
laboratory  angle  of  160  . 

The  targets  used  in  this  experiment  were  tritiated 

titanium  foils.   The  foils  were  5y  thick  and  the  activity 

2 
was  6  Curies/in  .   This  activity  corresponds  to  an  atomic 

ratio  of  approximately  one  between  the  tritium  and  the 
titanium.   The  tritium  contained  in  the  foils  was  deter- 
mined by  comparing  the  yield  obtained  by  the  solid  state 

3       3 
detector  with  the  known  yield  of   H(p,p)  H.   The  measured 

2 

tritium  values  for  the  two  targets  used  were  0.14  mg/cm 

and  0.11  mg/cm  .   The  error  in  these  concentrations  is 
determined  by  the  uncertainty  in  the  experimental  geometry 

and  is  estimated  to  be  10%.   The  errors  associated  with  the 

3       3 
statistics,  the   H(p,p)  H  cross  section  and  the  beam  charge 

measurement  combined  contribute  about  1%  to  this  error. 
Experimental  Data 
The  experimental  data  were  taken  in  four  series.   The 
first  series  consisted  of  a  set  of  excitation  curves  which 
were  utilized  to  determine  the  experimental  asymmetry  of  the 
angular  distribution.   The  second  series  consisted  of  angular 
distributions  measured  at  eight  angles  between  42   and  142  . 
These  measurements  were  extended  in  the  third  series  to  in- 
clude a  measurement  at  an  angle  of  20  .   The  fourth  series 
was  an  extension  of  the  excitation  curve  data  to  lower 
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energies  than  covered  in  the  first  series.   This  data  was 
then  analyzed  to  determine  the  effects  due  to  contributions 
from  the  electric  quadrupole  interaction. 

In  the  first  series,  the  excitation  curves  were 
measured  for  laboratory  angles  that  transform  approximately 
to  the  center  of  mass  angles  of  55°,  93°,  and  125°  (Figures 
2.5,  2.6).   The  measurement  was  taken  in  500  keV  steps  over 
the  energy  range  from  17-30  MeV.   The  resulting  cross  sections 
are  shown  in  Table  2.1.   The  center  of  mass  angle  in  the  for- 
ward direction  (9f  -  55°)  and  in  the  backward  direction 
(0B  ^  125°)  were  chosen  to  be  approximately  the  zeros  of  the 
second  order  Legendre  polynomial  P~  (cos  6).   This  condition 
on  6f  and  6B  requires  that  6B  =  180-6f  which  yields  the 
following  results  for  the  Legendre  polynomials: 

P1    (cos  6B)  =  -P1(cos  Gf  ) 

p„  (cos  eR)  =  p,(cos  ef  )  =  o 

Z  Z  (2.1) 

P   (cos  9B)  =  -P3(cos  6f  ) 

P4  (cos  9B)  =   P4(cos  6f  ) 

The  observed  yield  as  a  function  of  angle  can  be  written  as 

an  expansion  of  Legendre  polynomials, 

Y(8.)  =  AQ(1  +  £,akPk(cos  6)). 
k=l 


3       4 
Figure  2.5    H(p,y)  He  Center  of  Mass  Cross  Section  for 

9.  ,  =  90°.   The  data  from  this  series  of  experiments  is 
lab 

shown  as  the  solid  dots.   The  error  bars  represent  the 

statistical  error  associated  with  the  data  points.   The 

data  is  plotted  both  as  a  function  of  the  laboratory  energy 

of  the  incident  proton  and  the  center  of  mass  energy  of  the 

emitted  gamma  ray.   The  experimental  cross  section  reported 

13 

by  Meyerhof,  et  al.    is  included  as  the  open  squares.   The 

solid  curve  is  a  least  squares  fit  of  the  data  in  the  proton 

energy  range  from  5  to  30.5  MeV  to  a  cubic  polynomial.   The 

cross  section  was  obtained  by  normalizing  to  the  data  of 

13 
Meyerhof.      The  measured  target  thickness,  solid  angle 

and  collected  charge  was  combined  with  the  normalization  to 

obtain  the  detector  efficiency.   The  value  of  27  ±  3%  obtained 

is  consistent  with  the  efficiency  26  ±  6%  reported  by  Weller, 

4-   t  82 
et  al. 


3       4 
Figure  2.6     H(p,y)  He  Center  of  Mass  Cross  Sections  for 

6    =55   and  125°.   The  experimental  data  is  shown  as  the 

cm  c 

solid  dots  with  the  error  bars  indicating  the  associated 
statistical  error.   The  data  is  plotted  both  as  a  function 
of  the  laboratory  energy  of  the  incident  proton  and  the 
center  of  mass  energy  of  the  emitted  gamma  ray.   The  solid 
lines  are  fits  to  the  data  in  terms  of  the  model  discussed 
in  Chapter  III. 
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Under  the  assumption  that  multipolarities  of  L  >  2  contri- 
bute negligibly  to  the  cross  section,  the  sum  of  the  yields 
for  the  angles  8^  and  9   reduces  to 

Y(6f)  +  Y(6B)  =  2Ao(l+a4P4(cos  6f )  )  .  (2.2) 

where  a   and  a.  are  defined  in  equation  (1.4).   The  dif- 
ference between  the  yields  is  then 

Y(8  )-y(0B)  =  AQ(a1P1(cos  ef)+a3P3(cos  6f)-a1P1(cos  6g) 

-a3P3(cos  9B))  (2.3) 

=  2A  (a-,P1(cos  6f)+a3P3(cos  0f)). 

Therefore,  assuming  the  condition  a,P.(cos  8^)<<1,  which 
should  be  expected  if  the  E2  intensity  is  small  compared  to 
the  El  intensity,  the  total  cross  section  is  given  by 

a.     =    4ttA   c:  27t(Y(0.)+Y(6t,))  (2.4) 

u         O  LB 

and 

AQ  -  h    (Y(6f)+Y(6B))  .  (2.5) 

Under  these  conditions,  the  experimental  asymmetry  in  terms 
of  the  yields  Y(6f)  and  Y(6  )  is  equivalent  to  that  defined 
in  terms  of  the  coefficients  of  the  Legendre  expansion  of 
the  angular  distribution  (Figure  2.7).   This  equivalence  is 
shown  by  the  following  relations: 

Y(0f)  -  Y(0B)        =  2Ao(a1P1(cos  6f)+a3P3(cos  0f)) 


P1(cos  0f)  (Y(6f)+Y(0B))    P1(cos  0f)  (2Ao(l+a4P4(cos  6f ) ) 

( 2  fi ) 
b  a   +  P3(cos  0f)a3  t         b  (cos  e  )«1. 

p1(cos  ef) 


T         4  „, 

Figure  2.7     H(p,y)  He  Angular  Distributions.   The 
experimental  data  is  shown  as  the  solid  dots  with  the 
error  bars  indicating  the  associated  statistical  error. 
The  solid  lines  are  least  squares  fits  to  the  data  in 
terms  of  the  Legendre  Polynomials  as  discussed  in 
Chapter  III.   The  measured  values  for  the  cross  sections 
are  included  in  Appendix  B. 
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v/(e)o 


v/(e)o 
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As  will  be  shown  in  Chapter  III,  the  measured  asymmetry  and 
the  measured  cross  section  can  be  used  to  obtain  an  experi- 
mental E2  cross  section. 

The  El  cross  section  is  obtained  by  observing  that 


the  terms  of  the  multipole  expansion  (equation  1.5)  for 

da 
dil 


multipolarities  up  to  L  =  4  reduce  to  -^-r-  =  B  (a  +  1)  at 


0    =  90°. 
cm 

Therefore,  taking  a<<  1  (no  spin-flip  or  magnetic 
dipole  radiation)  the  yield  at  this  angle  is  the  El  cross 
section.   '   '     The  angle  measured  was  8,  .  =  9  0   so  that 

the  results  could  be  readily  compared  with  those  of  previous 

13,56 
experimenters . 

A  second  experimental  series  was  conducted  to  con- 
firm the  energy  behavior  of  the  E2  cross  section  and  to 
provide  a  direct  measurement  of  the  E2  amplitude.   This 
series  consisted  of  a  set  of  angular  distribution  measure- 
ments taken  every  2  MeV  between  17  and  27  MeV  incident 
proton  energy,  and  included  angular  distributions  measured 
for  proton  energies  of  13  and  30  MeV.   For  each  angular 
distribution,  yields  were  measured  for  a  minimum  of  seven 
angles  in  the  range  42  -142   covered  by  the  detector. 
Additional  angles  were  included  as  needed  to  define  the 
shape  of  the  peak  and  the  slope  of  the  sides  of  the  angular 
distribution.   The  minimum  number  of  measurement  points  was 
chosen  so  that  the  number  of  degrees  of  freedom  associated 
with-  a  Legendre  polynomial  least  squares  fit  to  four 
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coefficients  is  greater  than  zero.  The  amount  of  charge 
collected  was  varied  so  that  a  minimum  of  200  counts  was 
observed  for  the  142°  measurement.  The  yield  was  larger 
for  all  other  angles  observed  (Figure  2.7). 

The  third  series  of  measurements  was  conducted 
using  the  10.2  x  17.8  cm  Nal  detector.   The  purpose  of 
this  experimental  series  was  to  determine  the  shape  of  the 
angular  distribution  as  the  angle  approaches  zero.   Measure- 
ments were  made  every  2  MeV  in  the  proton  energy  range 

17  <  E   <   27  MeV.   Yields  were  measured  for  three  laboratory 
-   p  - 

angles  including  20° ,  9  0°  and  an  intermediate  angle  corres- 
ponding to  an  angle  at  each  energy  for  which  the  yield  was 
measured  in  the  second  series.   These  measurements  were 
normalized  to  the  data  of  the  second  series,  and  an  effective 
yield  for  the  25.4  x  25.4  cm  detector  at  the  laboratory  angle 
of  20°  was  obtained.   This  yield  was  incorporated  into  the 
angular  distributions  used  to  determine  the  Legendre 
expansion  coefficients. 

The  fourth  series  of  measurements  extended  the 
asymmetry  data  to  lower  energy.   The  TUNL  FN  tandem  was  used 
to  measure  the  fore-aft  asymmetry  every  1  MeV  in  the  proton 
energy  range  from  8  <  E    <  15  MeV.   This  information  was 
obtained  for  comparison  with  the  published  results  of 
earlier  experimenters. 


CHAPTER  III 
DATA  ANALYSIS 


Kinematics 

The  kinematics  of  the  reaction   H(p,y)  He  for 
proton  energies  above  10  MeV  are  subject  to  relativistic 
corrections.   Therefore,  the  kinematics  were  explicitly 

treated  according  to  special  relativity  following  the  treat- 

84 
ment  by  J.  J.  Jackson.     The  reaction  coordinates  in  the 

laboratory  and  center  of  momentum  frames  are  shown  in 

Figure  3.1. 


Laboratory 


Center  of  Momentum 


Figure  3.1   Laboratory  and  Center  of  Momentum  Coordinates 

The  conservation  equation  for  the  square  of  the  four  momen- 
tum of  the  alpha  particle  in  its  rest  frame  and  in  the 
laboratory  system  is 


/Z* 


(0,j   m  c  )   =  (  'y2   _!  m  ck-hk,  j   (y  m  c2+m.c2-E  ))2   (3.1) 
'  c  ^      P  /„PPty 
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where  y ■  is  the  Lorentz  transformation  parameter.   For  the 

c  

proton  in  the  laboratory,  this  parameter  is  Y   =     l-v2/c2 

■*        c  p  p 

Solving  the  conservation  equation  for  the  laboratory  gamma 
ray  energy  as  a  function  of  laboratory  angle  gives 

EY      5 ^ (3.2) 

2 (y  m  c2+m,c2-  ^y2-!  m  c2cos0  ) 
P  P     t      'p    p      y 

where  the  energy 

E"  =  (m^c't+m2c't  +  2Ypmpc2mtc2)J2  (3.3) 

is  the  total  energy  in  the  center  of  mass  system. 
The  gamma  ray  energy  in  the  center  of  mass  is 

(KM'^C  _  (34) 

'       2E' 

The  excitation  energy  of  the  compound  nucleus  is  found  by 
subtracting  the  alpha  particle  rest  energy  from  the  total 
center  of  mass  energy 

Ex  -  E'  -mac2  (3.5) 

The  relationship  between  the  angle  in  the  center  of 
momentum  and  in  the  laboratory  is  given  by 


i  _ 


1  Y  -1 

-  ,sin  6,       .   ,     ,.     'p         ,_  _, 

tan   (_  _j_ )  where  3  =   -^ —  ,     (3.6) 

Y   (cos6  -3)  Y2 

cm      Y  p 

where  the  Lorentz  transformation  parameter  for  the  center  of 

momentum  y    is  related  to  the  total  laboratory  enerqy  and 
cm  J      r,J 
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the  total  energy  of  the  compound  nucleus  by 

v    =  y  m  c2  +  m.c2  , -.    .,> 

' cm    'p  p t_  (3.7) 

E' 

These  relations  were  incorporated  in  the  computer  programs 
used  for  data  analysis. 

The  relationship  between  the  measured  cross  section 
in  the  laboratory  reference  frame  and  that  of  the  center  of 
momentum  reference  frame  is  also  necessary  for  comparison  of 
the  experiment  to  theory.   This  relation  can  be  obtained 
from  the  equality  of  the  yield  in  the  two  frames  following 
Marion.  This  equality  is 


^lab  <6Y)d"  =  Qcm  (e')dft1.  (3.8) 

Using  the  kinematic  relations  between  6  and  6 '  with  the 

definitions  of  the  solid  angle  increments,  which  are 

6.Q   =  2tt  sine  d9   and  dQ, '  =  2tt  sinG'dG'  ,  the  relation  be- 
Y   Y 

tween  the  two  differential  cross  sections  is  found  to  be 
as  follows: 

/nx    rr       /ni\  sin38'  ,cos8'cos9       .  , -,  Q. 

°lah{\)    =  acm(9  >  sTHe—  (  ...  .   T+  W"  (3'9) 

1  Y    smB  sm8 

1  Y 

This  form  of  the  relation  between  the  two  frames  was  con- 
venient to  use  as  the  analysis  programs  calculated  6 '  and 
Y    for  measured  gamma  ray  angle  9  . 
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Excitation  Curve  Analysis 

Assuming  only  El  and  E2  radiation  with  the 
additional  assumption  that  the  El  and  E2  spin-flip  terms 
are  negligible,  the  cross  section  may  be  written  in  terms 
of  the  spherical  harmonics  Y, ,  and  Y„, .   The  interaction  is 
modeled  by  assuming  that  it  consists  of  direct  and  resonance 
terms.   The  center  of  mass  cross  section  expression  is 


da    (6)        i6i  iS2      i($i+ni) 

— *jg =|ejDie    YXi  +  (  e2D2e    +  Rxe 


i(^2+Tl2) 

+R2e         )Y21 


(3.10) 

41 


This  expression  for  the  cross  section  implicitly  includes 
the  relations  a-.  =  -a-,  and  1  +  a„  +  a.  =  0  between  the  coef- 
ficients in  the  Legendre  expansion.   The  first  term  of  the 
cross  section  is  due  to  a  direct  El  interaction.   The 
second  term  is  due  to  a  direct  E2  interaction,   The  factors 
e,  and  e„  are  the  kinematic  effective  charges.     The  third 
and  fourth  terms  are  resolved  Breit  Wigner  resonances  cen- 
tered around  energies  E,  and  E*  respectively.   Each  term  of 
this  expression  will  be  examined  further  in  the  following 
paragraphs . 
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The  expressions  for  the  resonance  terms  in  the 
center  of  mass  are 


Rx  -  g^;)372  pVp  (3.ii) 


((£'  -  E  )2  +  r12/4)T 


and 


3/     1/ 


R2  =  g2(E')  '2  P  (3>12) 


((£'  -  E  )  2+r2/4)  2 
Y         2 

The  energy  dependence  of  the  resonance  partial  widths  is 
explicitly  included  in  these  expressions  where  P   is  the 
proton  penetrability  calculated  by  the  method  of  B.  Buck, 
et  al .       and 


2L+1        h 
3/   *  y'       k     v 

Y        k2  Ey 

describes  the  energy  dependence  of  the  transition  amplitude 
(T  )  of  a  L  =  2  multipole.   The  resonance  strength  for  the 
gamma  ray  resonant  absorption  is  then 

.■Vc'UL+l)   r/p  .  4^'(E;>2L+1PP  ,3.13, 

(E')!  (E')! 

The  energy  dependence  of  the  total  width  of  the  resonance 
is  included  by  using  the  proton  and  neutron  penetrabilities 
and  their  respective  reduced  widths,  i.e. 


r  =  r  +  r  =  2y2p  +  2y2p  .  (3.14) 

n     p     'n  n     'p  p 
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The  deuteron  and  gamma  partial  widths  are  assumed  to  be 
negligible  with  respect  to  the  neutron  and  proton  partial 
widths.   The  reduced  widths  are  solutions  to  the  equations 


r(E)  =  2v2P  (E)  +  2v2P  (E)  (3.15) 

'  n  n         PP 

and 

2   —    2 


for  the  resonance  centered  about  the  energy  E.  The  neutron 
and  proton  contributions  are  here  assumed  to  differ  only  by 
coulomb  effects.  The  phase  $  in  each  case  is  the  resonance 
phase,  given  by 

$.  =  tan"   (E  "  Ei)  (3.16) 

1  r±/2 

The  phases  n-,  and  ru  account  for  the  relative  phases  be- 
tween the  resonances  and  the  direct  terms. 

The  functional  form  for  the  direct  terms  was  taken 

Q  Q 

from  the  analysis  of  Flowers  and  Mandl.     They  assume 
gaussian  type  wave  functions  and  obtain  relations  of  the 
form 

±&  ,3/   -bE   i  6X 

D,  e     =  ai(E')V  e    P  e  (3.17) 


and 


16  3/o  s,      -bE    i62 

D2e      =  a2(E|)/2E  A  e    P  e     .  (3.18) 
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The  gaussian  wave  functions  parameters  are  chosen  to  fit 

the  4He  rras.  radius  of  1.61  fm.  given  by  electron  scatter- 

5/  k  -5 

ing  giving  a±   =    .101  (mb/MeV/2)^,  a2  =  5.13  x  10 

(mb/MeV1  ^2)ii   and  b  =  .043  MeV-1.41   For  comparison  to  the 
experimental  data,  the  direct  El  term  was  replaced  with  an 
empirical  function  and  the  strength  of  the  E2  term,  a2,  was 
allowed  to  vary.   The  phases  &±   and  62  were  varied  to  ob- 
tain the  relative  phase  between  the  direct  El  and  E2  terms 
that  best  fit  the  data. 

The  above  expression  for  the  differential  cross 
section  is  in  the  form  of  products  of  spherical  harmonics. 
These  products  are  expanded  in  terms  of  single  spherical 
harmonics.   The  product  expansion  is 


*"    *""  *.|»  -r|  M--L  '   4"<2L+1)    ]  (3-19) 

ooo   m'mM      ffl  ,, 

where  the  Clebsh-Gordon  coefficient  C™l£L   is  zero  for 
M  f   m'-m.   The  resulting  expansion  is  in  terms  of  the  YLQ 
functions  alone.   The  functions  YLQ  are  simply  related  to 
the  Legendre  polynomials  by  the  function 

YL0(e,*)  =  /^  PL  (cose)  (3.20) 

Collecting  the  terms  of  the  polynomials  and  comparing  with 
the  Legendre  expansion  of  the  angular  distribution  allows 
identification  of  the  expansion  coefficients  in  terms  of  the 
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direct  and  resonance  contributions  of  the  model.   In  order 
to  make  this  identification  the  following  three  definitions 
are  made: 

Bo  =  UjDj)2 

B      =    [2e1D1e2D2cos (6i-62)-2£iDi (Ricos ($r +n 1 -6 i ) +R2cos  ($2  +  n 2-6 i ) ) 

B      =    [  (e2D2)  2-2e2D2  (Ricos  ($!  +  ni-<S2)+R2Cos  ($2  +  ri2-62)  )     +Ri2 

+  2RxR2cos  (<J>  i+n  i  -<!>  2~n  2 )    +   R22]« 

(3.21) 

In  these  definitions  B   appears  in  terms  of  only  El  ampli- 
tudes, B-,  appears  in  terms  of  products  of  El  and  E2 
amplitudes,  and  B2  appears  in  terms  of  E2  amplitudes  only. 
The  coefficients  of  the  cross  section  expansion  may  then  be 
written  as  follows: 

Ao  =  va  <VB2)'  Aoai  =  T~z  Bl'  Aoa2  =  ^(_v!B2)' 

4  7T  v  D 


-12 
0~"      o  *  '   o      7(4tt;  (3.22) 


A„a3  -  -Aoai,  AQa,  =  n^y^2 


In  this  form  the  El  and  E2  contributions  are  explicitly 
included. 

The  Legendre  expansion  gives  the  differential  cross 
section  as  a  function  of  angle  6  in  the  limit  of  infinitesmal 
solid  angle;  however,  the  actual  detector  solid  angle  is 
finite.   In  order  to  compare  the  model  to  experiment,  it  is 
therefore  necessary  to  transform  to  finite  geometry.   This 


6  4 

transformation  introduces  geometric  correction  factors  Qi 
(Ferguson  89  )   into  the  expansion.   The  resulting  expression 
for  the  experimental  cross  section  is 


\pa(e)  =A  (Q  +aiQiPi (cos6)+a2Q2P2 (cos0) 

"M °   °  (3.23) 

+a3Q3P3  (cos6)+a4Q4Pi,  (cos8)  )  . 


The  yield  from  the  3H(p,y)4He  reaction  is  converted 
to  the  photonuclear  reaction   He(y,p)  H  by  the  principle  of 
detailed  balance.   The  relationship  between  the  cross 
sections  is  given  by 

a  (e)    1((2I(^H)+1)(2I(3H)+1))  (!V!!P}     (e).90 

Y,P       2      (2I(4He)+l)  E*      P'Y 

(3.24) 


where  1(A)  is  the  nuclear  spin  of  nucleus  A. 

Using  these  relations,  a  computer  program  was  de- 
veloped which  plotted  the  model  cross  sections  against  the 
measured  cross  sections  as  a  function  of  energy  and  angle. 
The  model  parameters  for  the  curves  that  best  represented 
the  plotted  experimental  data  were  used  to  obtain  the 
resonance  parameters  and  the  integrated  strength  of  the 
various  components  of  the  model.   These  results  are 
presented  at  the  end  of  this  Chapter. 
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Angular  Distributions 

The  angular  distribution  data  were  measured  in  an 
attempt  to  obtain  a  direct  measurement  of  the  E2  strength 
through  a  determination  of  the  coefficient  a^ .   The  magni- 
tude of  the  cross  section  for  6=0°,  180   also  provides  a 
determination  of  the  contribution  of  spin-flip  interactions 
to  the  cross  section.   The  experimental  data  are  fit  by 
linear  least  squares  to  the  first  four  terms  of  the  Legendre 
cross  section  expansion  for  finite  geometry  as  discussed  in 
the  previous  section.   For  some  measurements,  the  resulting 
fits  were  negative  for  6  =  0°,  180°.   A  non-linear  least 
squares  search  with  the  additional  requirement  that  the  cross 
section  be  positive  was  then  used  to  obtain  the  values  for 
the  coefficients  A  ,  a1#  a2,  a3#  and  a4 .   The  error  associa- 
ted with  each  coefficient  is  then  obtained  from  the 
curvature  of  the  Chi-square  space. 

The  assumption  that  the  El  and  E2  spin  flip  is 
negligible  is  tested  by  the  angular  distributions.   The 
truncated  Legendre  expansion  reduces  to 

g  (e)  =  A0(l+a1+a,+a,+a,)  (3.25) 

at  6  =  0°  and 

|2.  (9)  =  Ao(l-a1+a2-a3+alt) 

at  6  =  180°.  The  cross  section  is  equal  at  these  two  angles  for 
a3   =  -a:  and  zero  if  in  addition  l+a^a^  =  0.  Assuming  only  El 
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and  E2  Radiation,  the  coefficients  A  a,  and  A  a   can  be 

o  1       o  J 

written  in  terms  of  the  complex  vectors  |j  j R J>  for  total 

angular  momentum  (j-j)  coupling  where  j   =  i      +  S   is  the 

■+  -> 

total  angular  momentum  of  the  incident  particle,  j   =  S„ 

is  the  total  angular  momentum  of  the  target  and  J  =  ]   +  ]„ 
is  the  total  angular  momentum  of  the  system.   The  nuclear 
spins  are  S   and  SR  of  the  projectile  and  target  respec- 
tively/ while  I   is  the  orbital  angular  momentum.   The  R 
matrix  elements  (<I  X.   J  |  |  T  |  |  j  j  „J> )  for  radiation  of  multi- 
polarity  ££    and  compound  nucleus  of  ground  state  spin  1=0 
are  the  following  for  El  and  E2  radiation: 


i1^  i$3/ 

<011||T||W>  =  P^e    H,    <011|  |T|  I  3/21/2l>  =  P3/2  e   /2 


(3.27) 


l'$3,  1     95/ 

<022|  |T|  |  3/2  V22>    =    d3/2e  h  ,    <022|  |t|  |5/2V22>    =    d^e  /2  . 

The   resulting    expression   for    the   coefficients    are   as    follows: 

i(2$3/  -x*i,  )) 
Aoal   =   I   {/I°    Re(V3/2      e 

-/5  ,  ,  1{H%~H%\  (3.28) 

—  Re(P3/2d3/2      e  ) 

(2$5/  "^Sy   ) 


+  3/30    _     ,_       ,         „  /2  /; 


t-  Re(pv2d5/2  e  ]> 
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2  *     .        _    1  ,7, 


Aoa3   =   |   {/I5   Re(P^d3/2    e  /2  /2    ) 


+i^Re(P3/d3/    .*'    '*      "%    ) 

->  /2       /2 


(3.29) 


9^-n  i(1$3/-2*s/) 

— = —  Re(P3/  d5/  e  )  } 

->  n    /2 


The    relationship    a,    =    -a,    can    be    obtained    from    these    equa- 
tions  by    imposing    the    following    conditions    on    the    phases 


J$3/      =     l  $1/      +    nTT 

h  n 


n    =    0,2,4, 


(3.30) 


5/2 


3/ 

/2 


+  nit 


With  these  conditions  there  are  two  possible  classes  of 
solutions  to  equations  3.28  and  3.29  that  lead  to  the 
result  a.,  =  -a,  .   One  class  is  obtained  when  the  ratio  be- 
tween the  magnitudes  d3/   and  d5.   equals  -/6/3  in  which 

/2        h 

case  the  ratio  between  the  magnitudes  Pi.  and  P3   is 

/2  /2 

indeterminate.   The  other  possible  class  of  solutions  is 
obtained  when  the  ratio  between  the  magnitudes  P,.   and  P3 

/2  /2 

equals  -/2/2  in  which  case  the  ratio  of  the  magnitudes  d 3. 
and  d5/  is  indeterminate. 

/2 

Alternatively  the  system  may  be  described  in  terms 
of  the  (L-S)  coupling  vectors  ( I  £  S  £0S  LSJ>).   The  total 

Q  d  p  p 

anqular  momentum  is  t,   =  1      +  ~ta    and  the  total  spin  is 
r  a     p  ,       " 

S  =  £   +  §0.   The  total  angular  momentum  is  then  J  =  L  +  S, 
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These  vectors  can  be  related  to  the  ]  •  3  coupling  vectors 
of  equations  3.27  by  the  us. 
relationship  is  as  follows: 


91 
of  equations  3.27  by  the  use  of  the  9j  symbols.  '   The 


I   hlDh-LSJ>      =   J   [(2S+1)  (2L+1)  (2j  +1)  (2j  +l)]h 

a      p  ■  •  a       p 

JaJ3 


k        k        S  f 

(3.31) 


'     -e      [  I  V6J>' 

The  results  for  total  angular  momentum  and  parity  J   =  1 


are 

I  Oh    \h    101   >   =    -/3~    I    ¥sl    >   +    S5    I     V2V2I    > 
~T  3 

(3.32) 

1 0^  13s  ill   >  =     ft  I   hhl  >  +  /I  I    V2V2I   > 
~3  3 

IT  + 

The  results  for  total  angular  momentum  and  parity  J   =  2 


I  03s    2k    202    >    =    -/TO"         sA/22    >    +    ^    I      A/22    > 
"IT  5 

(3.33) 

I  Ois    23s    212    >    -       /15    I     V2V22    >    +    /H    I     5/2722       • 

The  two  cases  leading  to  the  relation  a3  =~a1  are  obtained 
by  either  |  0^1^111>  or  |  0^2J5212>  being  equal  to  zero. 

Studies  with  polarized  protons  indicate  that  |  0JslJslll>  is 

92  •  15,19,33 

non-zero.     As  seen  by  previous  experimenters 
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the  contribution  of  the  spin  flip  terms  |0J5l^lll>  and 
|.0^2Js212>  to  the  un-polarized  cross  section  is  zero  to 
within  the  experimental  error.   The  data  in  this  experiment 
are  consistent  with  neglecting  spin-flip  as  can  be  seen  in 
Figure  2.7  and  3.3.   Therefore,  these  terms  were  not 
included  in  the  analysis. 

Results 

The  experimental  data  is  consistent  with  the 

assumption  that  the  reaction  is  dominated  by  the  electric 

dipole  radiation.   Considering  only  electric  dipole  and 

quadrupole  radiation,  the  magnitude  of  the  cross  section 

for  0    =  90°  is  a  measure  of  the  electric  dipole  strength 
cm 

if  the  angular  distribution  of  the  quadrupole  radiation  is 
of  the  form  of  the  spherical  harmonic  Y„,  (which  vanishes 
at  90°).   Using  these  assumptions,  an  empirical  equation 
for  the  electric  dipole  strength  was  derived  from  the  ex- 
perimental data,  using  a  linear  least  squares  fit  to  the 
experimental   H(p,y)  He  cross  section  between  23  and  42  MeV 
excitation.   The  cross  section  in  this  energy  range  was 
found  to  be  best  represented  by  a  cubic  polynomial  with  a 
reduced  chi-square  of  0.02.   The  resulting  expression  for 
the  El  radiation  in  the  model  is 

_3  _4  _5   2  -7   3 

?    4m  c2E  {1.518x10   +3.953x10   E  -3.921x10   E  +5.419x10   E  } 
D   =    p    p  x  x  x 

E2 
Y  (3.34) 
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This  expression  gives  the  90   cross  section  shown  as  the 
solid  curve  in  Figure  2.5.   Using  this  result  for  D,  ,  the 
model  parameters  were  obtained  to  produce  the  solid  curves 
in  Figure  2.6  for  the  angles  0     =  53°,  123°.   The  asymmetry 
(a  )  resulting  from  this  calculation  is  displayed  as  the 
solid  curve  in  Figure  3.2.   The  parameters  of  the  calculation 
are  shown  in  Table  3.1  for  the  resonances  observed  in  the 
data. 

Table  3.1 
Direct  and  Resonance  Parameters  of  the  E2  Cross  Section 


E  (Mev)  * 

r     (keV) 

T     (MeV) 

IT 

n 

6 
i 

,mb         ,  ^ 

a( irr-) 

MeV  /2 

52 

35. 
39.5 

2.6 
0.60 

12.8 
3.2 

50° 
35° 

-65° 

1.24x10" 

-46° 

*all  energies  are  cm 

The  least  squares  fits  to  the  angular  distributions  are  shown 
in  Figure  2.7  as  the  solid  lines.   The  coefficients  for  the 
Legendre  polynomials  as  a  function  of  the  incident  proton 
energy  at  the  center  of  the  target  are  displayed  in  Table  3.2. 
The  Legendre  coefficients  obtained  from  the  excitation  curve 
analysis  are  shown  as  the  solid  curves  on  Figure  3.3  overlaid 
with  the  experimental  results.   The  experimental  results  are 
shown  to  support  the  hypothesis  that  a.,  =  -a,  and  l+a2+a4=0  to 
within  experimental  error  at  all  measured  energies  with  the 
exception  of  E   =23.0  MeV.   The  agreement  between 

IT 


Figure  3.2   Asymmetry  of  the  3H(p,Y)4He  Angular  Distribution. 
The  asymmetry  data  from  the  measurements  at  6cm  =  55   and 
125°  are  the  solid  dots  with  the  error  bars  showing  the 
statistical  error.   The  asymmetry  from  the  angular  Distri- 
bution of  Figure  2.7  are  shown  as  x's  with  arrowheads 
indicating  the  error  propagated  from  a.,  and  a3>   The  results 

i  4-+-  a    13,60,58 
reported  by  previous  experimenters  is  also  plotted. 
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Figure  3.3    Coefficients  of  the  Legendre  Expansion  of 
a(8).   The  a.  coefficients  from  the  angular  distributions 

of  Figure  2.7  are  shown  as  x's  in  the  upper  plot.   The 
error  bars  are  the  errors  of  the  coefficients  a.^  obtained 

from  the  least  squares  fit  of  the  data  to  the  Legendre 
expansion.   The  lower  plot  shows  the  experimental  AQ 

obtained  from  the  cross  sections  of  Figure  2.6.   The 
solid  curves  in  both  plots  are  the  results  obtained  from 
the  model  (Equation  3.10).   The  energy  E    is  in  the  center 
of  mass. 
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the  experimental  data  and  the  model  of  direct  E2  and  two  E2 
resonances  superimposed  on  a  dominant  El  background  inter- 
action is  seen  to  be  quite  satisfactory. 

For  comparison  to  theoretical  calculations  in  the 
literature,  the  integrated  cross  sections  due  to  the  indi- 
vidual components  of  the  model  were  evaluated.  The  total 
photonuclear  cross  section  is  obtained  by  integrating  the 
Legendre  polynomial  expansion  of  the  angular  distribution 
(Equation  1.4)  over  all  solid  angled,  The  result  is 


/dfl   da   =  4ttA  =  Bq  +  B2  =  a   +  a 


where  B   and  B~  are  defined  by  Equation  3.21.   As  seen  in 
Equation  3.21,  the  individual  contributions  of  the  direct 
and  resonance  terms  used  in  the  model  may  be  evaluated 
separately.   The  results  of  this  evaluation  are  shown  in 
Table  3.3. 
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Table  3.3 
Cross  Sections  Integrated  From  22  to  44  MeV  Excitation 


fa  dE  =  3  0.6    mb  MeV      /  aE1  dE  =  2  9.6  mb  MeV 

Y  /  P 


foE2      dE  =   1.01   mb  MeV     /  aE2  dE  =  8.12x10  "  mb/ 


MeV 


2 


E 

/   Ri2      dE   =      0.226    mb   MeV  /      Rl     dE    =    1.78xl0~      mb/ 

—  MeV 

E2 

2 


/  R2-   dE  =   0.04  7  mb  Mev      /    2  dE  =  3.07x10    mb/ 

—  MeV 

E2 

/ad    dE  =   0.969  mb  MeV      /a,    dE  =  6.69xl0~   mb/ 

d  MeV 


SUMMARY  AND  CONCLUSION 

The  TUNL  Cyclo-Graaff  was  used  to  obtain  data  from 
the   H(p,y)  He  reaction  for  protons  of  energies  from  8-30 
MeV  in  the  laboratory.   The  data  were  taken  as  a  series  of 
excitation  curves  and  angular  distributions.   The  data 
obtained  can  be  adequately  accounted  for  if  the  radiation 
is  viewed  as  being  predominantly  electric  dipole  with  an 
admixture  of  direct  and  resonant  electric  quadrupole  radia- 
tion.  Two  electric  quadrupole  resonances  are  observed  and 
the  resonance  width  obtained  for  each.   The  resonance  at  a 
compound  nucleus  excitation  of  35  MeV  with  a  width  of  12 
MeV  is  in  good  agreement  with  the  result  reported  by  Malcom, 
et  al.    in  the   He  photoneutron  experiments  at  the 
Saskatchewan  Linear  Accelerator  Laboratory.   A  second, 
previously  unobserved,  resonance  is  seen  with  a  width  of 
^3.2  MeV  at  an  excitation  energy  of  3  9.5  MeV.   A  survey  of 

the  experimental  literature  shows  a  possible  indication  of 

4        2 
structure  at  this  energy  in  the   He(y,d)  H  reaction  measure- 

69         4        3 
ments  reported  by  Skopik  and  Dodge    and  in   He(y/P)  H 

19 

assymetry  measurements  by  Arkatov,  et  al.,    thus  lending 

credibility  to  the  result  reported  here. 
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The  observed  resonances  are  in  good  agreement  with 

the  (J77,!1)  =  (2+,0)  and  (2+,l)  levels  reported  by  Szydlik, 

ofi  93 

et  al.    for  shell  model  calculations  with  modified  Sussex 

94 
potentials,  where  T  is  the  isospin  quantum  number.     This 

reference  quotes  levels  of  (2+,0)  located  at  about  34  MeV 
and  (2+,l)  located  at  about  44  MeV  excitation.   The  inte- 
grated strengths  of  these  two  levels  are  reported  as  0.86 
and  0.96  mb  MeV  respectively  giving  a  total  E2  strength  of 
1.86  mb  MeV  from  threshold  to  50  MeV.   This  strength  com- 
pares with  0.23  and  0.047  mb  MeV  for  the  two  levels  observed 
and  0.97  mb  MeV  for  the  E2  direct  term.   The  total  E2 
strength  from  this  experiment  is  found  to  be  1.01  mb  MeV 
(Table  3.3) . 

The  above  reference  would  imply  that  the  newly 
observed  level  at  39.5  MeV  would  have  an  isospin  assignment 
of  T  =  1.   A  definite  assignment  of  the  isospin  quantum 

number  can  not  be  made  from  this  experiment  as  both  T  =  0 

3       4 
and  T  =  1  quadrupole  states  can  be  formed  in  the   H(p,y)  He 

reaction.   The  narrow  width,  however,  suggests  the  T  =  1 

assignment  since  fewer  decay  channels  are  available. 

A  recent  microscopic  continuum  calculation  by 

95  + 

Bevelacqua  and  Philpott    indicates  a  (2  ,0)  resonance 

4 
of  3.8  MeV  width  at  about  39  MeV  in   He.   Bevelacqua  and 

Philpott  compare  this  level  to  the  experimental  data  of 

Lin,  et  al.96   Using  inelastic  scattering  of  deuterons  from 
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He  and  ^Li,  Lin  and  Associates  suggest  the  existence  of  a 
level  of  width  from  2.8  to  5.6  MeV  at  an  excitation  energy 
of  31.7  MeV  to  fit  the  trend  of  their  data  at  the  upper  end 
of  their  observed  energy  range.   Since  the  data  of  Lin,  et  al. 
includes  only  the  low  energy  tail  of  this  resonance,  their 
width  assignment  of  this  level  must  remain  questionable. 
While  it  is  tempting  to  associate  this  narrow  level  with  the 

level  observed  at  39.5  MeV  in  this  experiment,  the  level  at 

22 
35  MeV  observed  both  here  and  by  Malcom,  et  al.    cannot  be 

ignored.   If  the  level  in  the  reported  calculation  of 

Bevelacqua  and  Philpott  corresponds  to  the  35  MeV  resonance 

observed  here,  then  their  level  corresponding  to  the  39.5 

MeV  resonance  would  be  at  an  energy  beyond  the  range  of  their 

calculation. 

It  is  apparent  that  further  theoretical  work  is 

necessary  in  order  to  explain  and  understand  the  existence 

of  the  2+  level  at  39.5  MeV  observed  in  this  experiment. 

This  series  of  experiments  extends  the  energy  range  covered 

by  the  ^H(p,y)^He  reaction  measurements  from  18  to  30  MeV 

incident  proton  energy.   The  (p,Y)  reaction  lends  itself  to 

better  resolution  and  simpler  analysis  than  the  4He(Y/P)^H 

or  ^He(e,e'p)-^H  experiments  that  comprise  the  data  in  the 

literature  in  this  energy  range.   The  effects  of  an  E2 

component  in  the  radiation  are  seen  in  the  measured  angular 

distributions . 


APPENDIX  A 

The  spectra  for  each  energy  and  angle  measurement 
were  written  on  magnetic  tape  for  off  line  computer  analysis. 
The  spectra  were  fit  to  the  response  curve  of  the  detector 

for  monochromatic  y-ray.   This  response  curve  was  obtained 

3       4 
by  accumulating  a   H(p,y)  He  spectrum  with  2741  counts  in 

the  full  width  at  half  maximum  of  the  peak.   A  least  squares 

fit  to  the  spectrum  was  obtained  for  the  composite  function 

F(x)  =  (1-8  (x-C9  +  e)exp  (C  !+C2x+C  3x2+C\x3 ) 


+  9 (x-C9-e)exp  (C5+C6x+C7X2+C8x3) 


where 


r,  I  \  o   x  <  x^ 

6 (X-X   )   =   r  O 

ll   X  >  X 

-   o 

and  x   =  Cg  is  the  channel  where  the  two  functions  have  equal 
magnitude  and  slope.97  The  first  term  is  the  low  energy  tail 
of  the  peak  and  the  second  is  the  peak  and  high  energy  edge 
of  the  spectrum.   The  coefficients  C.  obtained  from  this  fit 
determined  the  standard  curve  shape.   The  experimental 
measurements  were  then  fit  to  the  equation 

F(X)  =B  (1-6 (X-C9+e)exp  (Ci+C2X+C3X2+C^X3) 
+  6(X-C9-e)exp  (C5+C6X+C7X2+C8X 3 ) 
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x-P 

where  X  = +  X    .   The  fitting  parameters  are  the 

R     max 

amplitude  B,  the  gain  R  and  the  peak  position  P.   The 
variable  x  corresponds  to  the  channel  numbers,  and  the  term 

X     is  the  peak  of  the  standard  spectrum.   The  fitting  is 

max        c 

done  by  a  parabolic  extrapolation  of  the  chi-squared  hyper- 
surface  routine  that  is  a  modification  of  the  routine  Chifit 

9  8 
described  by  Bevington.     This  routine  was  found  to  be 

reasonably  fast  while  yielding  reproducable  fits  to  the 

experimental  spectra.   The  coefficients  of  the  fitting 

equation  are  as  follows: 

3 

C,  =  3.6393  Cc  =  1.4269x10 

_2  1 

C,  =  8. 9808x10  C,  -  1.4674x10 

2  6 

1*  _2 

C   =  4.8943x10  C?  =  5.0242x10 

_7  _  5 

C4  =  9.1427x10  C8  =  5.7104x10 

X     =  312.5  CQ  =  294. 

max  9 

_  2 
R  =  (x2-Xi)  2.5x10         B   =  N/95. 

The  peak  channel  counts  N,  and  half  width  x2  -  x,  are 

initially  selected  by  the  operator  using  a  light  pen  on 

the  displayed  spectrum. 

A  sample  fit  is  shown  in  Figure  A.l  for  the  gamma 

ray  emitted  at  the  laboratory  angle  of  65   for  an  incident 

proton  energy  of  30  MeV. 


Figure  A.l    H(p,y)  He  Spectrum  for  30  MeV  Incident  Proton 
Energy.   The  full  energy  peak,  first  escape  peak,  and 
second  escape  peak  are  indicated.   The  solid  curve  indi- 
cates the  fit  obtained  using  the  peak  fitting  program 
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SiNDOD 


APPENDIX  B 

Table  B.l 
Experimental  Angular  Distributions* 


E  29.5           32.4           33.9           35.4 

Y 

6  26.  21.  21. 

o(6)  2.34  ±  .18     1.57  ±  .15     1.41  ±  .14 

9  44.            44.            44.            44. 

0(6)  4.95  ±  .20     4.48  ±  .15     4.42  ±  .15     4.14  ±  .14 

6  52.  52.  52. 

o(9)  5.63  ±  .17     5.93  ±  .17     4.80  ±  .16 

6  57. 

o(6)  7.59  ±  .25 

6  72.            73.            73.            73. 

0(6)  8.50  ±  .27     7.17  ±  .19     6.28  ±  .18     5.47  ±  .14 

6  82. 

o(6)  7.99  ±  .27 

6  92.            93.            93.            93. 

0(6)  6.70  ±  .24     6.02  +  .18     5.19  +  .17     4.78  ±  .27 

6  112.           113.           113.           113. 

o(6)  6.19  ±  .24     4.86  ±  .16     4.36  +  .16     3.85  ±  .15 

6  127.           132.           132.           127. 

o(8)  4.28  ±  .20     2.80  ±  .13     2.69  ±  .12     2.50  ±  .12 

6  143.           144.           144.           144. 

0(6)  2.17  +  .14     1.74  ±  .07     1.74  ±  .10     1.33  ±  .07 
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Table  B.l  (continued) 


E 
Y 

36.8 

38.3 

39.8 

42.0 

6 

21. 

21. 

21. 

0(6 

)              1.31 

+ 

.13 

1.42 

+ 

.14 

.86 

+ 

.11 

8 

44. 

44. 

44. 

0(6 

)              3.43 

+ 

.11 

3.68 

+ 

.11 

3.18 

+ 

.11 

6 

50. 

50. 

51. 

51. 

0(6 

)              3.86 

+ 

.12 

3.96 

+ 

.12 

3.45 

+ 

.11 

4.09 

+ 

.18 

8 

58. 

58. 

58. 

61. 

a(6 

1              4.80 

+ 

.13 

4.70 

+ 

.13 

3.92 

+ 

.12 

4.66 

+ 

.14 

0 

73. 

73. 

73. 

68. 

o(6 

6.68 

+ 

.17 

6.15 

+ 

.17 

5.00 

+ 

.15 

4.71 

+ 

.14 

0 

84. 

a(8 

5.33 

+ 

.22 

6 

93. 

93. 

93. 

97. 

0(0 

4.34 

+ 

.13 

3.42 

+ 

.11 

3.29 

+ 

.11 

3.38 

+ 

.19 

e 

113. 

113. 

113. 

113. 

o(6] 

3.24 

+ 

.11 

2.63 

+ 

.10 

2.86 

+ 

.11 

3.01 

+ 

.12 

0 

128. 

128. 

128. 

128. 

0(8] 

2.28 

+ 

.10 

1.72 

+ 

.08 

1.68 

+ 

.08 

2.01 

+ 

.10 

0 

144. 

144. 

144. 

144. 

0(0] 

1.32 

+ 

.06 

.99 

+ 

.04 

1.04 

+ 

.05 

.96 

+ 

.07 

all  entries  are  center  of  mass,  cross  sections  are  yb 
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